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lO ' Abstract 

We prove existence of a countable family of spherically symmetric self-similar wave 
\jj I maps from 3 + 1 Minkowski spacetime into the 3-sphere. These maps can be viewed 

\^ • as excitations of the ground state wave map found previously by Shatah. The first 

CS , excitation is particularly interesting in the context of the Cauchy problem since it 

plays the role of a critical solution sitting at the threshold of singularity formation. We 
analyze the linear stability of our wave maps and show that the number of unstable 
0^ . modes about a given map is equal to its excitation index. Finally, we formulate a 

^^ I condition under which these results can be generalized to higher dimensions. 

^ : 1 Introduction. 

Wave maps, defined as harmonic maps from a spacetime (M, r/) into a Riemannian manifold 
{N,g), have been intensively studied during the past decade (see the recent review [|l|). 
i^ ! The interest in wave maps (sometimes called also sigma models) stems from the fact that 

they contain many features of more complex relativistic field models but are simple enough 
to be tractable rigorously. In particular, the investigation of questions of global existence 
and formation of singularities for wave maps can give insight into the analogous, but much 
more difficult, problems in general relativity. With this motivation we have recently studied 
numerically the development of singularities for wave maps from 3 + 1 Minkowski spacetime 
into the 3-sphere P]. In this case it was known that: (i) solutions with small initial data 
exist globally in time p, |]; (ii) there exist smooth initial data which lead to blow-up in 
finite time. An example of (ii) is due to Shatah who constructed a spherically symmetric 
self-similar wave map of the form u{r,t) = fo{^^)- This solution is perfectly smooth for 



t < T but it breaks down at t = T. Our numerical simulations strongly suggest that 
the self-similar blow-up found by Shatah is generic in the sense that there is a large set 
of initial data which comprise the basin of attraction of the solution /q. In particular, 
it seems that all initial data of nonzero degree (which by definition are not small in the 
sense of [^ Q) blow up in this universal self-similar manner. The dynamical evolution 
of degree zero wave maps is more interesting because, depending on the "size" of initial 
data, the solutions either exist globally in time converging to the vacuum (this scenario is 
usually referred to as dispersion), or blow-up in finite time (where, as before, the blow-up 
profile is given by /o). Thus, in this case there arises a natural question of determining the 
boundary between the basins of attraction of these two generic asymptotic behaviors. In P| 
we studied this question numerically by evolving various one-parameter families of degree 
zero initial data interpolating between blow-up and dispersion. A typical initial data in 
this class is a gaussian with varying amplitude. We found that the initial data lying on 
the boundary between the basins of attraction of the solution /o and the vacuum solution 
converge asymptotically to a certain codimension-one attractor which is self-similar. This 
suggested that, besides /o, the model admits another self-similar solution, call it /i(j^), 
which has exactly one unstable direction. This expectation was confirmed numerically 
in [0. In a sense, the solution /i can be thought of as the excitation of the ground state 
solution /q. 

The aim of this paper is to give a rigorous proof of existence of a countable family of 
spherically symmetric self-similar wave maps from Minkowski spacetime into the 3-sphere. 
The above mentioned solutions /o and /i are the first two elements of this family. The 
proof is based on a shooting technique very similar to the one used by us in the case of 
harmonic maps between 3-spheres p. 

2 Preliminaries. 

A wave map f/ is a map from a spacetime M with metric t] into a Riemannian manifold 
A^ with metric g which is a critical point of the action 

The associated Euler-Lagrange equations 

U^U^ + T^c{U)daU^d''U^ = (2) 



constitute the system of seinilinear wave equations, where F's are the Christoffel symbols 
of the metric g. In this paper we consider the case where M = IR^^^, 3 + 1 dimensional 
Minkowski spacetime, and N = S^, the unit 3-sphere. In polar coordinates on IR^''^^ and 
S^ the respective metrics are 



-dt"^ + dr"^ + r'^duj'^, (3) 



and 



g = du^ + sin^ {u)dn^, (4) 

where du'^ and dQ'^ are the standard metrics on S'^, and u G [0, n]. We consider spherically 
symmetric maps of the form 

U{t,r,uj) = {u{t,r),n = iu). (5) 

Then the action (1) reduces to 

L = i I Lu', + ul + ^^^] r'dt dr du, (6) 

and the corresponding Euler-Lagrange equation is 

2 sm{2u) ^ .^. 

- Utt + Urr + -Ur 5 = 0. (7) 

This equation is invariant under dilations: if u{t,r) is a solution of equation (7), so is 
Ux{t,r) = u{Xt, Xr). It is thus natural to look for self-similar solutions of the form 

«(t,r) = /(^), (8) 

where T is a positive constant. As mentioned in the introduction such solutions are im- 
portant in the context of the Cauchy problem for equation (7) since they appear in the 
dynamical evolution as intermediate or final attractors. Substituting the ansatz (8) into 
(7) we obtain the ordinary differential equation 

2^. sin(2/) 
P P^(l-P^ 

where p = r/(T — t) and ' = d/dp. For t < T we have < p < oo. 

It is sufficient to consider equation (9) only inside the past light cone of the point (T, 0), 
i. e., for p G [0, 1]. This constitutes the two-point singular boundary value problem with 
the boundary conditions 

/(0) = and /(1) = |, (10) 



/" + -/'--S?^=0. (9) 



which are dictated by the requirement of smoothness at the endpoints. Once a solution 
of equation (9) satisfying the conditions (10) is constructed, it can be easily extended to 
p > 1 1^. Note that solutions of (9) and (10) are the critical points of the functional 



E[f] 



2 /■/2 



// 



2cos^(/) ^ 

l-p2 



dp, 



;ii) 



which, as was pointed out by Shatah and Tahvildar-Zadeh []^, can be interpreted as the 
energy for harmonic maps from the hyperbolic space H^ into the upper hemisphere of 
S^. Shatah |^ showed that E[f] is bounded from below over the if^-space of functions 
satisfying (10) and attains an infimum at a smooth function /q, the ground state solution 
of equation (9). Independently, Turok and Spergel [§ found this solution in closed form 



/o = 2arctan(p). 



(12) 



The central result of this paper is 



Theorem 1 There exist a countable family of smooth solutions fn of equation (9) satisfying 
the boundary conditions (10). The index n = 0, 1, 2, ... denotes the number of intersections 
of fn{p) with the line f = 7r/2 (the equator of S^) on p E [0, 1). 

Before proving this theorem in the next section, we present now some numerical results. 
As will be shown below the solutions satisfying /(I) = 7r/2 form a one-parameter family 
with asymptotics /(p) ~ 7r/2 + h{p — 1) near x = 1, while the solutions satisfying /(O) = 
form a one-parameter family with asymptotics /(p) ~ ap near p = 0. The solutions /„ are 
obtained by a standard shooting-to-a-fitting-point method, that is by integrating equation 
(9) away from the singular points p = and p = 1 in the opposite directions with some 
trial parameters a and b and then adjusting these parameters so that the solution joins 
smoothly at the fitting point. The discrete set of pairs (a„, 6„) generated in this way and 
the energies characterizing the solutions /„ are shown below for n < 4. 
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Figure 1: The ground state solution /o and the first three excitations generated numeri- 
cally. The solutions /„ with n > were first discovered numerically by Aminneborg and 
Bergstrom 0. 

3 Proof of Theorem 1 

To prepare the ground for the proof of Theorem 1 we first discuss some basic properties of 
solutions of equation (9). It is convenient to use new variables defined by 



1 



P 



cosh X ' 



and h{x) = f{p) 



71 



(13) 



The range of x is from a; = (for p = 1) to x = oo (for p = 0). Note that the number of 
intersections of / with the line / = tt/2 is the same as the number of zeros of h. In these 
new variables equation (9) becomes 



h" - coth{x)h' + sm{2h) = 0, 
and the boundary conditions (10) translate into 
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h{0) = and h{oo) = ±— , 



(14) 



(15) 



where the ± sign in the last expression, obviously allowed by the reflection symmetry 
h — ^ —h, is introduced for convenience. 



Lemma 1 For any b there exists a unique global solution hb{x) to equation (I4) such that 

h{x) ~ bx'^ (16) 

as X -^ 0. 
Proof: Defining v = h', let us rewrite equation (14) as the system of two integral equations 

v{x) = — sinh(x) / ds, h{x) = / v{s)ds. (17) 

Jo sinli s JO 

Following the standard procedure we solve (17) by iteration, setting 

;(-+i)(a;) = -sinh(a;) r sm(2/i(")(g)) ^^^ h^^+^\x) = f v^''\s)ds. (18) 



' ' /o sinhs JO 

with the starting values /i^^-* = bx"^ and f '•^^ = 2bx. It can easily be shown that the mapping 
{h^'^\x) ^ v^'^\x)) — ^ {h^'^^'^\ v'^'^^^\x)) defined by (18) is contractive for any finite x, hence 
the sequence {h^"'\v^'^^) converges to a solution of equation (14). The proof of uniqueness 
is also routine so we omit it. 

Definition 1 A solution of equation (I4) starting at x = with the asymptotic behavior 
(16) will be called the 6-orbit. Without loss of generality we assume that b > 0. The b-orbit 
which satisfies h{oo) = ±n/2 will be called a connecting orbit. 

Remark 1: In the following whenever we say "a solution" we always mean the 6-orbit. Also, 
when we say that some property holds for all x we always mean for all x > 0. We use lim 
to denote lima;^oo- 

Remark 2: The endpoints of connecting orbits {h = ±7r/2, h' = 0) are saddle-type critical 
points of the asymptotic {x —>■ 00) autonomous equation h" — h' + sm2h = 0. One can easily 
show (c/. IP) that the connecting orbits converge to these points along the one- dimensional 
stable manifolds ±h{x) ~ — 7r/2 + ae~^. 

The following function, defined for 6-orbits, will play a crucial role in our analysis 

W{x) = -h'^ + sin^h. (19) 

We have 

—- = coth(x)h'^, (20) 

dx 

so W is increasing (unless h is a. constant solution). Equations (19) and (20) imply that if 
W{xq) > 1 for some Xq ( and h is not identically equal to ±7r/2) then |/i'(2;)| > e > for 
X > Xq, hence limVr(x) = 00. Thus, if a 6-orbit crosses the line h = ±7r/2, then h' and h 
tend monotonically to ±00. 

6 



Lemma 2 A b-orbit (with nonzero b) which satisfies \h{x)\ < 7r/2 for allx, is a connecting 
orbit. 

Proof: We showed above that if W{xq) > 1 for some Xo, then \h\ tends to infinity, hence 
\h\ < tt/2 imphes that W{x) < 1 for all x, so liniW^(a;) exists. Thus, limW^' = which 
means by (20) that lim/i' = and next by (19) that limsin^/i exists, implying that also 
lim/i exists. By equation (14), sin2/i(oo) = since otherwise limh" ^ contradicting 
lim/i' = 0. Hence, h{oo) = ±7r/2 or h{oo) = 0. To conclude the proof note that the latter 
implies limW^ = which in view of (20) is possible only if W^ = 0, that is h = 0. 

The next two lemmas describe the behavior of 6-orbits for small and large values of the 
shooting parameter b, respectively. 

Lemma 3 If b is sufficiently small then the solution hi,{x) has arbitrarily many zeros. 

Proof: Define h{x) = hb{x)/b. The function h satisfies 

h" - cotHx)h' + ^^ = (21) 

with the asymptotic behavior h{x) ~ x^ as a; ^ 0. As 6 ^ 0, the solutions of equation 
(21) tend uniformly on compact intervals to the solution of the limiting equation 

H" - coth{x)H' + 2H = (22) 

with the asymptotic behavior H{x) ~ x^ as a; ^ 0. The solution H{x) can be found in 
closed form in terms of the hypergeometric function but for the purpose of the argument it 
is enough to observe that H{x) is oscillating at infinity, since this implies that the number 
of zeros of hb{x) = bh{x) increases to infinity as b tends to zero. 

Lemma 4 If b is sufficiently large then the solution hiy{x) increases monotonically to oo. 

Proof: As in the proof of Lemma 3, we use a scaling argument. This time, we define 
h{x) = hb{x/Vb). The function h satisfies 

h"-^coth(^]h'+'^ = (23) 

Vb \Vb) b 

with the asymptotic behavior h{x) ~ a;^ as a; ^ 0. As b ^ oo, the solutions of equation 
(23) tend uniformly on compact intervals to the solution of the limiting equation 

H" - -H' = 0, (24) 

X 



that is to H{x) = x^. Thus, on any compact interval the solution hi,{x) stays arbitrarily 
close to hx"^ if h is sufficiently large. In particular, hb{x) strictly increases up to some Xq 
where hi^xo) = 7r/2. Since W{xo) > 1, by the argument following (20) hf, tends monotoni- 
cally to oo. 

Now we are ready to prove Theorem 1. The proof will the immediate corollary of the 
following proposition 

Proposition 1 There exists a decreasing sequence of positive numbers {6„}, n = 0, 1, 2, ..., 
such that the corresponding bn- orbits are connecting orbits with exactly n zeros for a; > 0. 
Morever, lim„^oo ^n = 0. 

Proof: The proof is based on an inductive application of the standard shooting argument. 
Let Sq = {b\ hb strictly increases up to some xq where hb{xo) = 7r/2}. Let bo = inf 5*0. By 
Lemma 4 the set 5*0 is nonempty and by Lemma 3 bo > e > 0. The bo-OThit cannot cross 
the line h = 7i/2 at a finite x because the same would be true for nearby 6-orbits with 
b < bo, violating the definition of 6o- Thus, the 6o-orbit stays in the region \h\ < 7r/2 for 
all X, and therefore due to Lemma 2 it is a connecting orbit. By definition the 6o-orbit has 
no zeros for x > 0. 

To make the inductive step we need one more lemma. 

Lemma 5 If b = bQ — e for sufficiently small e > 0, then the solution hb{x) increases 
up to some xq where it attains a positive local maximum h{xo) < tt/2 and then decreases 
monotonically to — oo. 

Proof: By the definition of 69 there must exist a point Xq where /^^(aJo) = 0. Since by (14) 
a solution h cannot have a local minimum if /i > 0, it follows that there must be a point 
Xi > Xq where hh{xi) = (otherwise the 6-orbit would contradict Lemma 2). The idea of 
the proof is to show that W{xi) > 1 provided that e is sufficiently small. As argued above 
this implies that for x > Xi h^ decreases monotonically to —00. In the following we drop 
the index b on hb. From (19) we have 

^x^ phi xc\ ) 

W{xi) - W{xo) = / coth{x)h''^dx > - / h'dh. (25) 

Jxo Jo 

In order to estimate the last integral note that for x > Xq 

W{x) - W{xo) = -h'"^ + sin'^hix) - sm'^h{xo) > 0, (26) 



so —h' > Y2(sin h{xQ) — sin h). Inserting this into (24) gives 

W{xi) > - sin^hixo) + / J2{sm^h(xo) - sin^h) dh. (27) 

2 JO 

The right-hand side of this inequahty is an increasing function of h{xQ) which exceeds 1 
if 7r/3 < h{xQ) < 7r/2, as can be checked by direct calculation. The value hi,{xQ) will fall 
into that interval if e is sufficiently small because by continuous dependence of solutions 
on initial conditions, hb{xQ) -^ 7r/2 as e ^ 0. This concludes the proof of Lemma 5. 

Having Lemma 5 we return now to the proof of Proposition 1. Let 5*1 = {b\hb{x) 
increases up to some xo where it attains a positive local maximum h{xo) < vr/2 and then 
decreases monotonically up to some xi where h{xi) = —tt/2}. Let bi = inf 5*1. Due to 
Lemma 5 the set 5*1 is nonempty and by Lemma 3 6i is strictly positive. Using the same 
argument as above we conclude that the 6i-orbit must stay in the region \h\ < 7r/2 for all 
X, so it is a connecting orbit (asymptoting to —n/2). By definition the 6i-orbit has exactly 
one zero for x > 0. 

The subsequent connecting orbits are obtained by the repetition of the above construc- 
tion. Since the sequence {6„} is decreasing and bounded below by zero, it has a nonnegative 
limit. Suppose that lim„^oo^n = b* > 0. The 6*-orbit cannot leave the region \h\ < 7r/2 for 
a finite x because the set of such orbits is clearly open. Thus, the 6*-orbit is a connecting 
orbit with some finite number of zeros. But this contradicts the fact that the number 
of zeros of 6„-orbits increases with n. We conclude therefore that lim„^oo &« = 0. This 
completes the proof of Proposition 1. Returning to the original variables /(p) and p, and 
using the notation hn{x) = hb„{x), we have fn{p) = hn{x) + n/2 with /«(!) = 7r/2 and 
/n(0) = 0{modTT), as claimed in Theorem 1. 

We end this section with a remark about the large n limit. From lim„^oo bn = 0, it 
follows that lim„^oo ^n(3;) = for any finite x. The limiting solution h* = (or /* = 7r/2) 
is a singular map which geometrically corresponds to the map into the equator of S^. The 
"energy" of this map E[f*] = provides the upper bound for the "energies" of critical 
points of (11) (we write "energy" is quotation marks to emphasize that the functional (11) 
is not the true conserved energy associated with the action (6)). As follows from the proof 
of Lemma 3, the behavior of connecting orbits with large n (and consequently small 6„) can 
be approximated by the solution of equation (22), namely /i„(x) ~ bnH{x) on x G [0,x„) 
where Xn tends to infinity as n -^ oo. This fact can be used to prove some remarkable 
scaling properties of connecting orbits in the limit of large n. 

9 



For example one can show that (see the table in Section 2) 

E^ 27r 

lim -^ = e^. (28) 

For more detailed discussion of this issue we refer the reader to where the analogous 
behavior in the case of harmonic maps between spheres was derived. 

4 Stability 

The role of self-similar solutions /„ in the evolution depends crucially on their stability with 
respect to small perturbations. This problem was analysed by us in [§] by mixed analytic- 
numerical methods. In particular, we provided evidence towards the conjecture that the 
solution /o is asymptotically stable and, as such, has an open basin of attraction. To make 
the results obtained in [^ rigorous is a formidable task. In this section we discuss the 
first (easy) step in achieving this goal, namely we determine the character of the spectrum 
of the linearized operator around the solutions /„. A somewhat different but equivalent 
version of the linear stability analysis was presented in 0. 

We restrict attention to the interior of the past light cone of the point (T, 0) and define 



the new time coordinate s = — In J(T — ty — r^. In terms of s and p, Eq.(7) becomes 

e^"* / 9« X 2 sin(2M) , , 

- 71 ^ie-^^s)s + Upp + -Up - / ; = 0. 29 

Of course, this equation reduces to Eq.(7) if the solution is self-similar, that is s- independent. 
Following the standard procedure we seek solutions of (29) in the form u{s,p) = fn{p) + 
w{s,p). Neglecting the O(w^) terms we obtain a linear evolution equation for the pertur- 
bation w{s, p) 

''' [e-''ws)s + Wpp + -wp - ^^^^ w = 0. (30) 



(l_p2)2V- -^^^ ■ -PP ■ p-P P^il-P') 

Substituting w{s,p) = e^^~^^^'^v{p) into (30) we get the eigenvalue problem 



;i -^d_ f,d_\ ^ 2(1 - p') c os(2/„ 
p2 dp \ dp I p2 



Av = {l-X')v, where A = -^^^f^f p^f + ^^^ " ^ ^, ^^^"^ ■ (31) 



Note that the principal part of the operator A is the radial Laplacian on the hyperbolic 
space H^. We consider this problem in the space of functions which are square- integrable 

10 



on the interval [0, 1] with respect to the natural inner product on H^, that is 

veL'{[0,l] P dp). (32) 

(1-P ) 

In this function space A is self-adjoint hence its spectrum is real. Both endpoints are of the 
limit-point type. Near p = the leading behavior of solutions of (31) is v{p) ~ p" where 
«(« + 1) = 2, so admissible solutions behave as 

f (p) ~ p as p ^ 0. (33) 

Near p = 1 the leading behavior is v{p) ~ (1 — p)^ where /3 = (1 ± vA^)/2 so eigenfunctions 
must have A^ > and behave as (up to a normalization constant) 

t;(p) ~ (l-p)5(i+l^l) as p^l. (34) 

All A^ < belong to the continuous spectrum of A. The case A = will be treated 
separately below. Note that this eigenvalue problem has the symmetry A — * —A (that is 
why we wrote A + 1 rather than A in the ansatz for w) . Each eigenvalue A^ > gives rise to 
an unstable mode which grows exponentially as e^''^'"*"^-''*. To find the eigenvalues we need 
to solve (31) on the interval p G [0, 1] with the boundary conditions (33) and (34). In 
we did this numerically (for ra < 4) by shooting the solutions from both ends and matching 
the logarithmic derivatives at a midpoint. For example, for n = 1 we got A^ ~ 28.448; 
for n = 2 we got \\ ~ 28.132, Ag ~ 3372.12. Our numerics strongly suggested that the 
point spectrum of the operator A around the solution /„ has exactly n positive eigenvalues 
A^ > (fc = 1, . . . , ra). Now, we will give a simple proof of this property. 

The proof is based on the observation that the solution with A = corresponds to the 
gauge mode which is due to the freedom of choosing the blowup time T. To see this, consider 



a solution fn{r/{T' — t)). In terms of the similarity variables s = — In w (T — t)^ — r^ and 

p = r/(T — t), we have 

In other words, each self-similar solution fn{p) generates the orbit of solutions of (29) 
parametrized by e. It is easy to verify that the generator of this orbit 

11 



solves (30), thus v^^Xg^{p) = p^\ - p^f'^ip) satisfies (31) with A = OQ. Since vf^^^Xp) has 
exactly n zeros on p G (0, 1) (because /„ has n extrema), it follows by the standard result 
from Sturm-Liouville theory that there are exactly n positive eigenvalues, as conjectured 
in 10. To summarize, we showed that the self-similar solution /„ has exactly n unstable 
modes, which means in particular that the solution /q is linearly stable. 

Remark. If we view the solutions /„ as harmonic maps from the hyperboloid H^ into S*^, 
then the eigenvalue problem (31) determines the spectrum of the Hessian of the energy 
functional (11) 

s:2 77iri' i/„. „.\ f I „2„./2 , ^ '^os(^z/„ 




S''E[U]iv,v) = -J^ A"+ ^jr ^^^Mp. (37) 



Within this approach the argument given above can be rephrased in terms of the Morse 
index. In particular, it implies that the Morse index of the solution /q is zero, in agreement 
with Shatah's result that /o is a local minimum of the energy functional (11). Note that 
in this context the gauge mode acquires a geometrical interpretation as the perturbation 
induced by K, the conformal Killing vector field on H^, 



v):2.„. = £Kfn, where K = py 1 — p^ d/dp . (38) 



gauge 



Aside, we remark that by solving (31) one can show that the Morse index of the singular 
map /* = 7r/2 is infinite. This fact could be probably used to give an alternative proof 
of Theorem 1 via Morse theory methods using the ideas of Corlette and Wald developed 
recently in the case of harmonic maps between spheres fT^ . 



5 Generalization to higher dimensions 

The proof of Theorem 1 is rather robust which suggests that the result can be generalized 
in various directions. One possibility, which will not be pursued here, is to consider more 
general nonconvex targets^. Another possibility is to consider the analogous problem in 
higher dimensions, that is wave maps U : M —* N, where M = IR"^"^^, m + 1 dimensional 

^We emphasize that A = is not an eigenvalue because Vgwuge{p) is not square-integrable at p = 1. 
Hovewer, this solution is distinguished from the rest of the continuous spectrum by the fact that it is 
subdominant at p — 1 (such a solution is sometimes referred to as a pseudo-eigenfunction) . 

^For example, one can easily verify that the proof of Theorem 1 goes through if the metric (4) is replaced 
hy g — du^ + s^{u)d^^ where the function s(u) satisfies the following conditions (cf. |iri|): (i) s(0) — and 
s'(0) = 1; (ii) s{u) is monotone increasing from w = up to some u* > where it attains a maximum. 
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Minkowski spacetime, and A^ = S*™, the unit m-sphere. At the same time one can relax 
the equivariance ansatz (5) by admitting the maps of the form 

Uit,r,u) = iuit,r),n = xiuj)), (39) 

where x is a homogeneous harmonic polynomial of degree / > 0. The ansatz (5) is the 
special case of (39) with / = 1. Assuming self-similarity we obtain the analogue of Eq.(9) 



pff , / ^- 1 , (m-3)p \ _ fesin(2/) 
p 1 — p^ j p2(l — p2 



r+(^^^ + ^^^^^T/'- 477^ = 0' (40) 



where k = l{l + m — 2)/2. As before, we want to construct smooth solutions on the interval 
[0, 1] satisfying the boundary conditions (10). Standard analysis of the behavior of such 
solutions at the endpoints yields that /(p) ~ ap^ near p = and /(p) ~ 7r/2+6(l— p)*^™""^)/^. 
Note that the latter implies that the desired smooth solutions can exist only if the dimension 
m is odd. Although Eq.(40) looks more complex than (9), the same change of variables as 
in (13) transforms (40) into 

h" - (m - 2) cot\i{x)h' + k sin(2/i) = 0, (41) 

which has the same form as (14) apart from the change of constant coefficients. Now, let 
us see which steps of the shooting argument from Section 3 are effected by this change of 
coefficients. Lemma 1 holds with the asymptotic behavior near a; = replaced by h{x) ~ 
^^m-i Lg]2imas 2 and 4 remain valid because their proofs are dimension independent. The 
only fact which is dimension sensitive is Lemma 3, because in m dimensions the limiting 
equation analogous to (22) reads 

H" - (m - 2) coi\i{x)H' + 2kH = 0, (42) 

so Lemma 3 is true iff the solution H{x) is oscillating at infinity, that is, 4k > {m — 2)^. 
This imposes the condition 

l>y^{m-2). (43) 

Under this condition the proofs of Lemma 5 and Proposition 1 remain basically unchanged, 
thus we have 

Theorem 2 For each oddm > 3 and I satisfying the condition (43), there exist a countable 
family of smooth solutions fn of equation (40) satisfying the boundary conditions (10). The 
index n = 0, 1, 2, ... denotes the number of intersections of fn{p) with the line f = n/2 on 
pG [0,1). 

13 



This theorem extends the recent result of Cazenave, Shatah, and Tahvildar-Zadeh [^ who 
proved existence of the ground state solution /q in odd dimensions under the condition (43) 
using variational methods. 
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